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Abstract 

Continuing a line of investigation initiated by Texier and Zumbrun 
on dynamics of viscous shock and detonation waves, we show that 
a hnearly unstable Lax-type viscous shock solution of a semilinear 
strictly parabolic system of conservation laws possesses a translation- 
invariant center stable manifold within which it is nonlinearly orbitally 
stable with respect to small n perturbatoins, converging time- 
asymptotically to a translate of the unperturbed wave. That is, for 
a shock with p unstable eigenvalues, we establish conditional stability 
on a codimension-p manifold of initial data, with sharp rates of decay 
in all . For p = 0, we recover the result of unconditional stability 
obtained by Howard, Mascia, and Zumbrun. 
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1 Introduction 

In this paper, we continue a line of investigation opened in [TZl, TZ2, TZ3, 
TZ4, SS, BeSZ] going beyond simple stability analysis to study nontrivial 
dynamics, and associated physical phenomena, of perturbed viscous shock 
waves in the presence of linear instability. The above-mentioned references 
concern Hopf bifurcation to time-periodic behavior associated with transition 
to linear instability arising through the passage from stable to unstable half- 
plane of a complex conjugate pair of eigenvalues of the linearized operator 
about the wave. See also [ZH, Zl] for discussion of (nonstandard, due to 
embedding in essential spectrum) bifurcations associated with the passage of 
a simple eigenvalue through zero. 

In the present work, in the spirit of [GJLS, Li] and other works outside 
the shock wave context, we consider the situation of a viscous shock sub- 
stantially after the onset of instability, i.e., with one or more strictly unstable 
but no neutrally unstable eigenvalues, and seek to describe the nearby phase 
portrait in terms of invariant manifolds and behavior therein. Specifically, 
for shock waves of systems of conservation laws with artificial viscosity, we 
construct a center stable manifold and show that the shock is conditionally 
(nonlinearly) stable with respect to this codimension p set of initial data, 
where p is the number of unstable eigenvalues. As discussed for example in 
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[AMPZ, GZ], such conditionally stable shock waves can play an important 
role in asymptotic behavior as metastable states. 

Consider a viscous shock solution u{x,t) = u{x), lim^^-too = u±, 
without loss of generality stationary, of a semilinear parabolic system of 
conservation laws 

(1-1) Ut + f{u)^. = u^^, 

/ e ]R", a; , i e ]R, under the basic assumptions: 
(HO) / e k>2. 

(HI) A± :— df{u±) have simple, real, nonzero eigenvalues. 
Linearizing (1.1) about u yields linearized equations 

(1.2) ut = Lu := -{df{u)u)^ - u^^, 

for which the generator L possesses [He, S, ZH] both a translational zero- 
eigenvalue and essential spectrum tangent at zero to the imaginary axis. 

The absence of a spectral gap between neutral (i.e., zero real part) and 
stable (negative real part) spectra of L prevents the usual ODE-type decom- 
position of the flow near u into invariant stable, center, and unstable mani- 
folds. The first result of this paper, by now little more than a remark, is that 
we can still determine center stable and unstable manifolds, and that these 
may be chosen to respect the underlying translation- invariance of (1.1). See 
[TZl] for closely related results on existence of translational- invariant center 
manifolds. As the needed ingredients do not seem to be found in one place, 
we nonetheless for completeness carry out the proof in full detail. 

Theorem 1.1. Under assumptions (HO)-(Hl), there exists in an neigh- 
borhood of the set of translates of u a codimension-p translation invariant 
(with respect to H^) center stable manifold M.cs, tangent at u to the cen- 
ter stable subspace of L, that is (locally) invariant under the forward 
time- evolution of (3.1) and contains all solutions that remain bounded and 
sufficiently close to a translate of u in forward time, where p is the (neces- 
sarily finite) number of unstable, i.e., positive real part, eigenvalues of L. 

Next, specializing a bit further, we add to (HO)-(Hl) the additional hy- 
pothesis that uhe a. Lax-type shock. 

(H2) The dimensions of the unstable subspace of df{u-) and the stable 
subspace of df{u^) sum to n -|- 1. 
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We assume further the foUowing spectral genericity conditions. 
(Dl) L has no nonzero imaginary eigenvalues. 

(D2) The orbit u{-) is a transversal connection of the associated standing 
wave equation Ux = f{u) — f{u-). 

(D3) The associated inviscid shock is hyperbolically stable, i.e., 

(1.3) det(rf , . . . , r~_i, r++i, . . . , r+, {u+ - U-)) ^ 0, 

where rf,...r~_i denote eigenvectors of df{u-) associated with negative 
eigenvalues and . . . r+ denote eigenvectors of df{u^) associated with 
positive eigenvalues. 

As discussed in [ZH, MaZl], (D2)-(D3) correspond in the absence of a spec- 
tral gap to a generalized notion of simplicity of the embedded eigenvalue 

A = of L. Thus, (D1)-(D3) together correspond to the assumption that 
there are no additional (usual or generalized) eigenvalues on the imaginary 
axis other than the transational eigenvalue at A = 0; that is, the shock is not 
in transition between different degrees of stability, but has stability properties 
that are insensitive to small variations in parameters. 

With these assumptions, we obtain our second and main result character- 
izing stability properties of u. In the case p = 0, this reduces to the nonlinear 
orbital stability resuh estabhshcd in [ZH, MaZl, MaZ2, MaZ3, Z2, Z3, Z4]. 

Theorem 1.2. Under (H0)-(H2) and (D1)-(D3), u is nonlinearly orhitally 
stable under sufficiently small perturbations in L} fl H'^ lying on the codi- 
mension p center stable manifold /Acs of u and its translates, where p is the 
number of unstable eigenvalues of L, in the sense that, for some «(•), 

\u{x,t) - u{x - a{t))\LP < C{1 +t)~^^^~p^\u{x,0) - u{x)\Linm, 
^^^^ \u{x,t) - u{x - a{t))\Hi < C{l + t)''^u{x,0) -u{x)\Linm, 

a(t) < C(l + t)-^\u(x, 0)-u(x)\LinH^: 
a{t) < C\u{x, 0) - u{x)\Linm- 

Moreover, it is orbitally unstable with respect to small perturbations not 
lying in A4, in the sense that the corresponding solution leaves a fixed-radius 
neighborhood of the set of translates of u in finite time. 
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1.1 Discussion and open problems 



It is easily checked that the results of this paper go through for general semi- 
linear parabolic term Buxx, B constant, under the standard assumptions of 
[ZH]. An interesting open problem is to extend these results to the quasi- 
linear and or partially parabolic ("real viscosity") case. We note that the 
only issue here is to establish existence of the center stable manifold, as the 
proof of stability goes through essentially unchanged, incorporating the nec- 
essary modifications detailed in [Z2, Z3] to deal with the quasilinear/partially 
parabolic case. See also related discussion in [Li]. 

Another interesting problem would be to extend our conditional stability 
result to the case of nonclassical under- or overcompressive shocks using 
pointwise estimates as in [HZ, RZ]; see Remark 4.7. 

Finally, we mention the problem of determining conditional stability of 
a planar standing shock u{x, t) = u{xi) of a multidimensional system of 
conservation laws 



which likewise (by the multidimensional arguments of [Zl, Z2, Z3]) reduces 
to construction of a center stable manifold, in this case involving an infinite- 
dimensional unstable subspace corresponding to essential spectra of the lin- 
earized operator L about the wave. 

Plan of the paper. In Section 2 we give for completeness a particularly 
concise proof of the center stable manifold theorem for ODE. In Section 3, 
loosely following [TZl], we show how to extend this to semilinear parabolic 
PDE, while preserving the key property of translation-invariance. Finally, in 
Section 4, we estabhsh conditional stabihty by a modification of the argu- 
ments of [Z4, MaZ2, MaZ3] in the stable (p = 0) case. 



The center stable manifold construction in the PDE case follows closely the 
construction for finite-dimensional ODE, which we therefore recall here for 
completeness; see also [B, VI]. Consider an ODE 
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2 Center Stable Manifold for ODE 



(2.1) 



u = 
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and an equilibrium /(w*) = 0, with associated linearized equation 
(2.2) v' = Av, A := df{u^). 

Associated with A, define the center stable subspace T^cs as the direct sum 
of all eigenspaces of A associated to neutral or stable eigenvalues, i.e., eigen- 
values with zero or positive real part. Likewise, define the unstable subspace 
S„ as the direct sum of eigenspaces associated to unstable eigenvalues, i.e., 
eigenvalues with strictly positive real part, so that M" = S^s © S^. 

Defining the associated (total) eigenprojections lies and as the sum of 
all eigenprojections associated with neutral-stable and unstable eigenvalues, 
respectively, we have, either by reduction to Jordan form or direct estimation 
using the inverse Laplace transform/resolvent estimates, bounds 

|e^*ne.| < C(r7)e^* t > 0, 

■ |e^*n„| < C(^)e-''l*l, t < 0, 

for any rj > strictly smaller than the minimum of the real parts of unstable 
eigenvalues and ^ > arbitrarily small. 

Proposition 2.1 (Center Stable Manifold Theorem for ODE). For f in 
(jk+i^ 1 < A; < OO; there exists local to it* a center stable manifold M.cs, 
tangent at to T^cs, expressible in coordinates w :— u — as a C'^ graph 
$cs • '^cs ^ Scs © that is (locally) invariant under the flow of (2.1) 
and contains all solutions that remain bounded and sufficiently close to in 
forward time. In general it is not unique. 

2.1 Frechet differentiability of substitution operators 

Standard invariant (e.g., stable, center, center stable) manifold constructions 
proceed by fixed point iteration on various weighted L°° spaces 

(2.4) 11/11, := supe^(*-*°)|/(t)|, B, {/ : ||/||, < oo}, 

t>to 

with T] positive for stable manifolds and negative for center or center stable 
manifolds. As described in [B, VI], Prechet differentiabihty of the associated 
fixed-point mapping (hence eventual smoothness of the resulting manifold) 
hinges on Frechet differentiability with respect to spaces of the special 
class of substitution operators, defined for g : ^ M" as 

(2.5) G{fm:=g{f{t)). 
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Let denote the Banach space of C'' functions : R" — > with \d^g\ 
uniformly bounded for < j < /c, with associated norm 

WaWc^ ■= Y] supld-^^l- 

0<j<k 

Lemma 2.2. For r] > 0. k > 1, if g e C^{W W), then G is from 
Br, ^ Br,, with d^Gf{t) = {d'g){f{t)). 

Proof. See Appendix A. □ 

Lemma 2.3. Let g e C^+^ and < —rj' < —r]/k. Then, the substitution 
operator G is G^ from Br,' Br,, with d^Gf{t) = {d^g){f{t)). 

Proof. More generally, the result of Lemma 2.3 holds for any g e C^"*"", 
< a < 1, in the sense that \d^g{x + h) — d''g{x)\ < G\h\°' for some uniform 
C > 0, as may be seen by rewriting the {k — l)th-order Taylor remainder 

formula g{x + h) - %-ig{x, h) ^ d^g{x + dh)^-^0^de^h^ as 

where Tkg{x, h) denotes the fcth-order Taylor expansion of g about x evalu- 
ated at x + h. then using the assumed uniform bound on \d^g[x + h) —d^g{x)\ 
to obtain \g{x + h) — Tkg{x, h)\ < C\h\^'^"', and thus 

\g{x + h)- Tkg{x, h)\^^k+aW < C|/i|J,+". 

Without loss of generality taking a sufficiently small, this yields the result. 
A similar estimate yields continuity of the kth Prechet derivative. □ 



2.2 Smooth dependence of fixed-point solutions 

We next present two general results on smooth dependence of fixed point 
solutions. Let T{x, y) be continuous in x, y and contractive in y, T : Bi x 
B2 B2 for Banach spaces Bi and B2, defining a fixed point map y{x), 
y : Bi ^ B2, continuous in the parameter x, such that y{x) — T{x,y{x)). 
Then, we have the following standard result. 
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Lemma 2.4. IfT is Lipschitz in {x. y). then y is Lip schitz from Bi -^32- 
IfT is (Frechet sense), k > 1, in {x, y), then y is C*^ frovn Ei B2, with 

(2.6) {dy/dx){xo) = (Id - Ty)-'T,{xo, y{xo)) 

and higher derivatives {d^y/dx^){xo), 1 < j < k obtained by formal differen- 
tiation of (2.6); substituting for lower derivatives wherever they appear. 

Proof. See Appendix A. □ 

The next Lemma shows how we can recover dependence of fixed point 
solutions in the case, as in Lemma 2.3, that T is differentiable only from a 
stronger to a weaker space. We discuss higher derivatives later where they 
appear, since they involve specific chains of successively weaker spaces that 
arc not convenient for statement as a general theorem. Let B'2 C B2, with 
II ■ ||b^ > II ■ ||b2) ^^cl T{x,y) be a map Bi x B2 ^ B'2 that is Lipschitz 
continuous in {x,y) and contractive (with respect to || • ||b^) in y. Denote 
by y{x) : Bi ^ B'2 the unique Lipschitz fixed-point solution defined by 
y{x) = T{x,y{x)). 

Lemma 2.5. If (i) T is continuously differentiable from Bi x B'2 ^ B2, and 
(a) Ty extends to a bounded linear operator from B2 — * B2, continuous in 
operator norm with respect to {x,y), with \Ty\1S2 < 1, then y is continously 
differentiable from Bi — > B2, with {dy/dx){xo) = (Id — Ty)~^Tx{xo, y{xo)) . 

Proof. By Taylor's Theorem and |||/(a;2) — y{xi)\\s'^ < L\\x2 — a;i||B^, we have 

y{x2) - y{xi) = T^{x2 - xi) + Ty{y2 - yi) 

(2.7) + o(||,T2 - xiWb^ + \\y{x2) - y{xi)\\B'^) 

= T^{x2 - xi) + Ty{y2 - yi) + o(||a;2 - Xi\\bJ, 

where the o{\\x2~Xi H^J term is measured in the weaker || ■ \\1s2 norm. Observ- 
ing by \Ty\js2^B2 < 1 s-iid Neumann series inversion that (Id — T^) considered 
as an operator from B2 B2 is invertible with uniformly bounded inverse 

|(Id - Ty)~^\B2^B2 < (1 - \Ty\B2^B2)~^ , 

we may solve (2.7) to obtain 

y{x2) - y{xi) = (Id - Ty)~^T^{x2 - xi) + o{\\x2 - a^iHsJ, 

yielding differentiability as claimed, with [dy/dx) — {ld — Ty)~^Ty. continuous 
by the assumed continuity of and Ty as operators from Bi — > Bi and 
B2 ^ B2. □ 
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2.3 Global Center Stable Manifold construction 



We now establish a global version of 2.1 for small Lipschitz nonlinearity. As 
in Section 2.2, denote by the space of functions that are uniformly 

bounded in up to k derivatives and consider for a fixed, constant matrix A 
and an arbitrary nonlinearity N such that N{t, 0) = 0, Nyj{t, 0) = the ODE 

(2.8) w' ^ Aw + N{t,w). 

Proposition 2.6. For N E C^^^ with Lipschitz constant s > sufficiently 
small, (2.8) has a unique invariant manifold A4cs tangent at w = to the 
center stable subspace 'Ecs of A, consisting of the union of all orbits whose 
solutions grow at sufficiently slow exponential rate \w{t)\ < Ce^'*' in positive 
time, for any fixed 9 < 9 < rj. 

Proof. Applying projections Uj, j — cs,u to (2.8), we obtain using the vari- 
ation of constants formula equations 

Ujw{t) = e^(*-*°'^)n^-'u;(toj) + / e^^^-'^UjN{s, w{s)) ds, 

Jto,j 

j = cs,u, so long as the solution w exists, with to.j arbitrary. Assuming 
growth of at most \w{t)\ < Ce^* in positive time, we find using (2.3) and 
the bound |iV(u')| < e\w\ coming from iV(0) = and the assumed Lipschitz 
bound on N, that as to,u — ^ +oo, the first term e^^*~*°'^^n„w(toj) converges to 
zero while the second, integral term converges to J^°° e^^*~^''IluN{s, w{s)) ds, 
so that, choosing to.cs = 0, we have 

Ucsw{t) = e^*ne.«;(0) + / e^(*-^)ne.A^(s, w{s)) ds, 

Jo 

/ + 00 
e'^^*-''>IluN{s,w{s))ds. 

Summing, we obtain for Wcs '■= llcs'w{0) the fixed-point representation 
w{t) = T{wc„ w) := e^'wcs + f e^^'-'^U,sN{s, w{s)) ds 
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valid for solutions growing at rate at most \w{t)\ < Ce^* for t >0. 
Define now the negatively-weighted sup norm 



.,-:=supe-''*|/WI, 



noting that \f{t)\ < e^*||/||_g; for all t > 0. We first show that, for \wcs\ and 
5,e > sufficiently small, the integral operator T is Lipschitz in Wcs and 
contractive in w from the || • ||_g[-ball B{0, 6) to itself. 

Using (2.3), \N{t,w{t))\ < e\w{t)\, and \w{t)\ < e^l*l||^w||_^, we obtain 

\T{w){t)\ < Ce^\'\\wcs\+Ce\\w\\_Q(^j\'\'-'^e^^'Us + £'^ 

hence, using 77 ± 6^ > and taking Cie < 1/2 and C\wo\ < 5/2, that 

\\T{w){t)\\_§ < C\wcs\ + Ce\\w\\_§(^ y* e^l*-^le^'(l^l-W) ds 

/+00 _ 
evit-s)^ei\s\-\t\) ^^j < ^i^^^i ^ Cie\\w\\_0 < S. 

Similarly, we find that 

\\T{wi)-T{w2)\\_g < Ce||^«i -«;2||_e~(y*e^l*-^le^>lds + jf^°°e''(*-^)e^>lds) 

< Cie\\Wi - W2\\_g < (1/2)||W1 - W2\\_§, 

yielding contraction on B(0, 5) and thus existence of a unique fixed point 
w = w(wcs)- A similar estimate shows that T is Lipschitz in Wcs, so that w{-) 
is Lipshitz from to B_q. 

We next investigate smoothness of w{-). Note that T decomposes into 
the sum of a bounded, hence C°°, linear map Wcs — > e'^^Wcs from Ecs — > B_g 
and the composition /C • jV' of a bounded (hence C°°) linear map 

/•t f+00 
W) I e^^'-^^U^sfis) ds - jf e^(*-)n„/(s) ds 

from B_Q B_Q and a substitution operator N'{w){s) := N{w{s)) with 
e C^+^ that by Lemma 2.3 is from B_q. Moreover, the 
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first derivative dM{w){s) = dN{w{s)), by \dN\ < s, extends as a bounded 
linear operator from B^qi B^0>, any 6' > 0, with \dN\is_g, < e, whence 
extends as a bounded Unear operator from B-e> B^g', that for any given 
9' > 0, in particular 9, is contractive, \Tyu\B_g, < 1, for £ > sufficiently 
small, independent of w and Wcs- Applying Lemma 2.5, we find that w{wcs) 
is from — > B^o' for any 9' > and £ > sufficiently small, with 

(2.10) d^/;-(Id-T^)-X... 

Differentiating (2.10) using the chain, product, and inverse linear opera- 
tor derivative formulae, validated by the fact that (Id — Tyj)~^ is uniformly 
bounded (since |T^| < 7 < 1) from B^gi — > B^0> for all < ^0 < ^' < ^70 < 
and the observation that higher (mixed) partial derivatives of T exist and are 
continuous from B-o'/^k+i) ^-0' on the same range, we find that w{wcs) 
is from B_g^^ B_g for C > sufficiently large, so long as 6*0 > is 

chosen < 9/C and e > is taken sufficiently small. 
Finally, defining 

r+00 

(2.11) ^wcs) := Uuw(ws)\t=o = - / e^(*-^)n„7V(s, w{s)) ds, 

Jo 

we obtain as claimed a function from J^cs ^u, whose graph is the invari- 
ant manifold of orbits growing at exponential rate \w{t)\ < Ce^* in forward 
time. Prom the latter characterization, we obtain evidently invariance in 
forward and backward time. By uniqueness of fixed point solutions, we have 
w{0) = and thus $cs(0) = 0. Finally, differentiating (2.11) with respect to 
Ws, we obtain by Lemma 2.3 

r+00 

d$(0) = - / e^^^-'^U^N^(s, 0)(dw/dws)(s) ds = 0, 
Jo 

since A^u,(0) = 0, yielding tangency as claimed. □ 



2.4 Local construction: proof of Theorem 2.1 

We can reduce the general situation, locally, to the case described in the 
global result by the following truncation procedure. Consider a general non- 
linearity N{t,w). Introducing a cutoff function 

p{x) = 
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define N%t,w) := p{\w\/e)N{t,w). 

Lemma 2.7. Let N e C''+\ k > 1, and N{t, 0) = 0, <9^iV(t, 0) = 0. Then, 
jye g C^'^^, = N for \u\ < e, and the Lipschitz constant for with 
respect to w is uniformly bounded by Ce for all < e < eo, where 

(2.12) C = 2(1 + max|r|p'(|r|)) max \dlN{t,w)\. 

Proof. See Appendix A □ 

Proof of Proposition 2.1. Defining w :— u — u^, we obtain the nonlinear per- 
turbation equation 

(2.13) w' = Aw + N{w), 

with A constant and N E C^'^^ satisfying N{0) = 0, dN{0) = 0. Applying 
the truncation procedure, we obtain a modified equation 

(2.14) w' ^ Aw + N%w) 

for which A^^ e C^^^ with arbitrarily small Lipschitz norm £ > and N = 
within a neighborhood B(0,e) oi w — 0, i.e., with identical local fiow. Ap- 
plying Proposition 2.6, we obtain a global center stable manifold for (2.14), 
which is therefore a local center stable manifold for (2.13). Noting that solu- 
tions that stay uniformly bounded and close to the equilibrium for positive 
time are also bounded solutions of the truncated equations, we find that all 
such belong to the constructed center stable manifold. □ 

Remark 2.8. The inclusion of t-dependence of N in (2.8) is not needed for 
the present application (2.13), but allows also the treatment of time-periodic 
solutions after Floquet transformation to constant-coefficient linear part. 

3 Center Stable Manifold for PDE 

We now turn to the case of a general semilinear parabolic equation 

(3.1) Ut — J-{u) :— h{u, Ux) + Uxx, u,h & R", 
with steady state u{x,t) = u{x) and associated group invariance 

(3.2) ^^{u){x,t) -.^uix + a^t), 
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and construct a local translation- invariant center stable manifold in the vicin- 
ity of u. 

Our approach follows closely that used to construct translation-invariant 
center manifolds in [TZl], by introduction of a reduced flow on the quotient 
space induced by group equivalence. However, we coordinatize differently, 
using orthogonal projection rather than eigenprojections, to avoid the diffi- 
culty (as in the ultimate application to shock waves) that the zero eigenvale 
associated with translation-invariancc may be embedded in essential spec- 
trum of the linearized operator about the wave, as a consequence of which 
there may not exist a well-defined zero eigenprojection with respect to 
(see [ZH] for further discussion) . 

We make the following assumptions, in practice typically satisfied [TZl]. 

(AO) h G C^+\ k>2. 

(Al) The linearized operator L = ^{u) about u has p unstable (positive 
real part) eigenvalues, with the rest of its spectrum of nonpositive real part. 
(A2) \diu{x)\ < Ce-^l^l, ^ > 0, for 1 < j < A; + 2. 

Proposition 3.1 (Center Stable Manifold Thm. for PDE). Under assump- 
tions (A0)-(A2), there exists in an neighborhood of the set of translates 
of u a translation invariant ( with respect to H^) center stable manifold 
M-cs, tangent at u to the center stable subspace T,cs of L, that is (locally) in- 
variant under the forward time-evolution of (3.1) and contains all solutions 
that remain bounded and sufficiently close to a translate ofu in forward time. 
In general it is not unique. 



3.1 Reduced equations 

Differentiating with respect to a the relation = dt{^°'{u)) — T{^°'{u)), we 
recover the standard fact that 

^ d^a{u) du 
da \oc=o dx 

is an zero eigenfunction of L, by the assumed decay of u^. 
Define orthogonal projections 

(3.3) n^r^^^, ni:=id-n2. 
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onto the range of right zcro-cigenfunction (f) :— [d/dx)u of L and its orthog- 
onal complement (j)^ in L^, where (•, •) denotes standard inner product. 

Lemma 3.2. Under the assumed regularity h e C'''^^, k > 2, Uj, j — 1,2 
are bounded as operators from to itself for < s < A; + 2. 

Proof. Immediate, by the assumed decay of = Mj. and derivatives. □ 

Introducing the shifted perturbation variable 

(3.4) v{x,t) :— u{x -\- a{t),t) — u{x) 

similarly as in [Z4, MaZ2, TZl], we obtain the nonlinear perturbation equa- 
tion 

(3.5) dtV^Lv + g{v)-dta{(t) + d^v), 
where L := ^{u) and 

(3.6) Q{v) = g{v,Vx,x) := h{u + v,Ux + v^) - h{u,Ux) - dh{u,Ux){v,Vx) 

is a quadratic-order Taylor remainder. 

Choosing dta so as to cancel 112 of the righthand side of (3.5), we obtain 
finally the reduced equations 

(3.7) dtv = U,{Lv + g{v)) 
and 

7r2(Lv + g{v)) 



(3.8) dta 



1 + 7r2(9^^;) 



for V G 0"*", where 7T2V := (0, f)||0|^2, of the same regularity as the original 
equations. 

Clearly, (3.8) is well-defined so long as \dxv\L°° < C\v\h2 remains small, 
hence we may solve the v equation independently of a, determining a- 
behavior afterward to determine the full solution 

u{x, t) — u{x — a{t)) + v{x — a{t),t). 
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Moreover, it is easily seen (by the block-triangular structure of L with re- 
spect to this decomposition) that the linear part IliL = riiLIli of the v- 
equation possesses all spectrum of L apart from the zero eigenvalue asso- 
ciated with eigenfunction 0. Thus, we have effectively projected out this 
zero-eigenfunction, and with it the group symmetry of translation. 

We may therefore construct the center stable manifold for the reduced 
equation (3.7), automatically obtaining translation-invariance when we ex- 
tend to the full evolution using (3.8). See [TZl] for further discussion. 

3.2 Preliminary estimates 

For ease of notation, introduce Lq :— HiL, Qq :— UiQ. 

Lemma 3.3. Under the assumed regularity h e C'^'^^, both Q and Qq are 

Frechet dijferentiable of order {k + 1) considered respectively as functions 
from to and (j)-^ C to H^: Q on the whole space and Qo for \v\h^ 
sufficiently small. 

Proof. Differentiability of Q follows by direct calcTilation; sec [S, TZl]. Dif- 
ferentiability of IIi^ follows similarly, using also the fact, already discussed, 
that 1 -|- 7T2{dxv) remains bounded from zero for \v\h2 small, and the fact (see 
Lemma 3.2) that 11^ as bounded linear operators from each to itself are 
infinitely differentiable in the Prechet sense. □ 

Lemma 3.4. Lq generates an analytic semigroup e^°* = Ilie''"*!!! on 0-*- C 
H^. Moreover, the unstable (positive real part) spectra of L and Lq agree in 
both location and multiplicity, with associated total unstable eigenprojections 
n° and Uu related by 11° = nin„ni and total center stable eigenprojections 
and lies related by U^^ = nillcslli. Likewise, except possibly at X = 0, 
the resolvent sets of L and Lq agree, with (A — Lq)~^ = Ili{X — L)~^Ili. 

Proof. Direct computation using relations LIli = L and LII2 = yields 
the resolvent relation, whence we obtain the remaining relations by their 
characterizations in terms of the resolvent (for example, the characterization 
of eigenprojection as residue of the resolvent operator [Kat]). As L is a 
sectorial operator, it follows that Lq is as well, and both generate analytic 
semigroups given by the inverse Laplace transform of the resolvent. □ 
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Corollary 3.5 ([TZl]). Under assumptions (A0)-(A2), 



for some P > 0, and for all cu > 0, for all t >0. 

Proof. These follow from the corresponding estimates for L, which are stan- 
dard semigroup estimates for second-order elliptic operators; see Appendix 
A. □ 

Let p be a smooth truncation function as in 

§2.4ande?o'(^) p{^)Oo{v). 

Lemma 3.6 ([TZl]). The map : x R ^ x R ts C^+^ and its 
Lipschitz norm with respect to v is 0{S) as S 0. 

Proof. See Appendix A. □ 
Corollary 3.7. Under assumptions (A0)-(A2), 

for some (3 > 0, and for all lu > 0, for all t >0, with Lipshitz bounds 

We'^^^n^^dg'.WH^^H^ < cj{i + rV2)e-*, 
We-'^-^n^dg^oWH^^H^ < cje-^K 

3.3 Translation-invariant center stable manifold 

Proof of Proposition 3.1. Using bounds (3.9), (3.10), and (3.11), and observ- 
ing that, since finite-dimensional, the unstable flow e^^^llu is well-defined in 
both forward and backward time, wc find, applying the finitc-dimcnsional ar- 
gument word-for-word, that the center stable manifold of the truncated fiow 
Ut = LoM-f^o('") graph of the map $ defined on H^s as ^{ucs) = (t'cs(O)), 
where Vcs is the unique solution in 0-*- C of 

pt poo 

v{x, t) = e*^°iie. + / e(*-^)-^°ne.e?o^(T;(x, s))ds - / e(*-^)-^°n„^?o^(T;(x, s))ds 
Jo Jt 

as guaranteed by the Contraction Mapping Theorem. 
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Indeed, the only deviation from the finite-dimensional case is the appear- 
ance of the new factor (1 + (t — s)^^^"^) in the integrand of estimates having to 
do with integral term fl^ e^^''^^~^^IlcsGo, which, since integrable, does not alter 
the final estimates. The proof of smoothness relies on these same bounds, so 
likewise carries over word-for-word as in the finite-dimensional case. Finally, 
invariance in forward time follows by the characterization of the center sta- 
ble manifold as the set of sohitions of the truncated equations growing no 
faster than \v{t)\ < Ce'^* in forward time. (Since the flow of (3.1) is only a 
semigroup, we cannot conclude invariance in backward time as in the finite- 
dimensional case.) Since sufficiently small bounded solutions of the original 
are bounded solutions also of the truncated equations, they are contained in 
the center stable manifold, independent of the choice of truncation function. 

This gives a center stable manifold with the stated properties for the 
reduced equation (3.7) for v E (p-^. Solving for the shift a in terms of v using 
(3.8), and substituting v, a into (3.4) to obtain u, we obtain a translation- 
invariant center stable manifold for the original equations (3.1). Observing 
that small bounded solutions v of the reduced equations, by (3.4), correspond 
to solutions u of the (3. 1) remaining close to a translate of u, we are done. □ 

3.4 Application to viscous shock waves 

Proof of Theorem 1.1. By Proposition 3.1, it is sufficient to verify that (HO)- 
(Hl) imply (A0)-(A2), for h{u,Ux) := f{u)^ = df{u)Ux. Clearly, (HO) implies 
(AO) by the form of h. Plugging u = u{x) into (1.1), we obtain the standing- 
wave ODE f{u)x — UxX, or, integrating from — oo to x, the first-order system 

Ux = f{u) - f{u-). 

Linearizing about the assumed critical points u± yields linearized systems 
Wt = df{u±)w, from which we see that u±. are nondegenerate rest points by 
(HI), as a consequence of which (A2) follows by standard ODE theory [Co]. 

Finally, linearizing PDE (3.1) about the constant solutions u = u^, we 
obtain Wt = L±w := —df{u±)wx — Wxx- By Fourier transform, the limiting 
operators L± have spectra \f{k) — —ikaf{k) — A;^, where the Fourier wave- 
number k runs over all of M; in particular, L± have spectra of nonpositive 
real part. By a standard result of Henry [He], the essential spectrum of L 
lies to the left of the rightmost boundary of the spectra of L±, hence we 
may conclude that the essential spectrum of L is entirely nonpositive. As 
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the spectra of L to the right of the essential spectrum by sectoriahty of L, 
consists of finitely many discrete eigenvalues, this means that the spectra 
of L with positive real part consists of p unstable eigenvalues, for some p, 
verifying (Al). □ 

4 Conditional stability analysis 

Define similarly as in Section 3.1 the perturbation variable 

(4.1) v{x,t) := u{x + a{t),t) - u{x) 

for u a solution of (1.1), where a is to be specified later in a way appropriate 
for the task at hand. Subtracting the equations for u{x + a{t),t) and u{x), 
we obtain the nonlinear perturbation equation 

(4.2) vt-Lv^ Niv%, 

where L := —dxdf{u) + 9^ as in (1.2) denotes the linearized operator about 
u and 

(4.3) N{v) := -U{u + v)- f{u) - df{u)v) 
where, so long as [t'li^i (hence \v\loc and ImI^oo) remains bounded. 



2^ 



(4.4) N{v) = 0{\v 
4.1 Projector bounds 

Let n„ denote the eigenprojection of L onto its unstable subspace E^, and 
Hcs = Id — the eigenprojection onto its center stable subspace Scs. 

Lemma 4.1. Assuming (HO)-(Hl), 

(4.5) n,a, = djlj 

for j — u, cs and, for all 1 < p < oo, < r < 4, 



(4.6) 



|ncs|iy.p->M/'''Pi |ncs|iy.p-»w""'P ^ c, 

\^cs\wr,p^W^-P , \^cs\wr,p^W^-P ^ C. 
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Proof. Recalling (see the proof of Theorem 1.1) that L has at most finitely 
many unstable eigenvalues, we find that may be expressed as 

p 

i=i 

where (f)j, j = l,...p are generalized right eigenfunctions of L associated 
with unstable eigenvalues A^, satisfying the generalized eigenvalue equation 
(L — XjY^4>j = 0, Vj > 1, and (pj are generalized left eigenfunctions. Noting 
that L is divergence form, and that Xj ^ 0, we may integrate {L — \jY^(f)j — 
over M to obtain A^^ J (j)jdx — and thus J (pjdx — 0. Noting that (f)j, and 
derivatives decay exponentially by standard theory [He, ZH, MaZl], we find 
that 

with $j and derivatives exponentially decaying, hence 
Estimating 

for 1/p+l/q — 1 and similarly for d^Iluf, we obtain the claimed bounds on 
Hu and Ilu, from which the bounds on lies = Id — Hu and lies = Id — !!„ 
follow immediately. □ 

4.2 Linear estimates 

Let Gcs{x, t; y) :— Ilcse^^Sy{x) denote the Green kernel of the linearized solu- 
tion operator on the center stable subspace Sc^. Then, we have the following 
detailed pointwise bounds established in [TZ2, MaZl]. 

Proposition 4.2 ([TZ2, MaZl]). Assuming (H0)-(H2), (Dl)-D(3), the cen- 
ter stable Green function may he decomposed as Gcs — E + G, where 

(4.7) E(x,t;y) ^d^u(x)ej{y,t), 
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(4.8) e{y,t) = J] (err/n (^^) - errfn (^:^)) l.iv) 
for y < and symmetrically for y > 0, /j^ e M" constant, and 

n 

\G{x,t;y)\ < Ce-''(l^-^l+*) + r^/2g-(x-y-a-t)VMtg-^x+ 

fe=i 

(4.9) MKt\>\y\r ^ ^ 

a^>0,aj<0 

I V f-1/2 -{x-a+{t-\y/a^\))^/Mt -^x- 

+ 2^ ^{|afc*l>l2/|}^ e ^ e 



aj, >0,aj>0 



(4.10) 

n 

A:=l 

I V ^-1/2 -(x-a-(t-|y/a^|))VMt -r;a;+ 



aj. >0, a. <0 



I V V ^-1/2 -(a;-at(t-|y/a^|))VMt 

+ 2^ ^{|afc*l>|j/ir e ^ '= e j 



aj. >0,aT>0 



for y < and symmetrically for y > 0, for some rj, C, M > 0, where 
are the eigenvalues ofdf{u±), denotes the positive/negative part ofx, and 
indicator function X{\a^t\>\y\} -'- /'^'^ l^fe^l — Ivl (^i^d otherwise. 

Proof. As observed in [TZ2] , it is equivalent to establish decomposition 

(4.11) G^Gu + E + G 

for the full Green function G{x,t;y) :— e^*Sy{x), where 

p 

G^{x,t;y) := n„e^*5,(x) = e^* ^ 0,(x)4(y)* 

for some constant matrix M e C^^'^ denotes the Green kernel of the linearized 
solution operator on S„, 0^ and right and left generalized eigenfunctions 
associated with unstable eigenvalues Xj, j — 1, . . . ,p. 
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The problem of describing the full Green function has been treated in 
[ZH, MaZ3] , starting with the Inverse Laplace Transform representation 

(4.12) G{x, t; y) = e^'dy{x) = j e^\\ - L{e))-^5y{x)d\ , 

where 

r:= d{X : < rji - r]2\^X\} 

is an appropriate sectorial contour, 771, 772 > 0; estimating the resolvent kernel 
G\{x, y) :— (A — L{e))~^5y{x) using Taylor expansion in A, asymptotic ODE 
techniques in x, y, and judicious decomposition into various scattering, ex- 
cited, and residual modes; then, finally, estimating the contribution of various 
modes to (4.12) by Riemann saddlepoint (Stationary Phase) method, moving 
contour F to a optimal, "minimax" positions for each mode, depending on 
the values of (x, y, t). 

In the present case, we may first move F to a contour F' enclosing (to the 
left) all spectra of L except for the p unstable eigenvalues Aj, j — 1, . . . ,p, 
to obtain 

G{x,t;y)= i e^*(A - L)-^(iA + VRcsiducA.(,)(e^*(A - L)-^5j,(a;)), 

where Residue ^^(e) (e^* (A — L)^^5y{x)) = Gu{x,t;y), then estimate the re- 
maining term f^, e^*(A — L)^^dX on minimax contours as just described. See 
the proof of Proposition 7.1, [MaZ3], for a detailed discussion of minimax es- 
timates E + G and of Proposition 7.7, [MaZ3] for a complementary discussion 
of residues incurred at eigenvalues in {JJA > 0} \ {0}. See also [TZl]. □ 

Corollary 4.3 ([MaZl]). Assuming (H0)-(H2), (D1)-(D3), 

/+00 
G{;t;y)fiy)dy\L. < G{1 + t)-'^^^-'^^\f\L., 
-00 

/+00 
Gy{;t;y)fiy)dy\Lr < C(l + i)-^^^-^)-^ 
-00 

for all t > 0, some C > 0, for any I < q < p (equivalently, 1 < r < p) and 
f G L'^, where 1/r + 1/q = 1 + 1/p. 
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Proof. Standard convolution inequalities together with bounds (4.9)-(4.10); 
see [MaZl, MaZ2, MaZ3, Z2] for further details. □ 

Corollary 4.4 ([Z4]). The kernel e satisfies 

\ey{;t)\L.,\et{;t)\L.<Ct-'^^'-'/^\ 

for all t > 0. Moreover, for y <0 we have the pointwise bounds 
\ey{y,t)\,\et{y,t)\ <Ct ^ 2^ (e Mt +e m* j, 

«fc>0 

|etj,(|/,t)| < Cr^ (^e M^+e s^j, 

a->0 

/or A'/ > sufficiently large, and symmetrically for y > 0. 

Proof. Direct computation using with definition (4.8); see Appendix A. □ 

4.3 Reduced equations II 

Recalling that dxU is a stationary solution of the linearized equations Ut — Lu, 
so that Ldx=0, or 

/oo 
G{x,t;y)ux{y)dy = e^^u^ix) = dxu{x), 
-oo 

we have, applying Duhamel's principle to (4.2), 

/oo 
G{x,t;y)vo{y) dy 
-oo 

pt POO 

— / Gy{x,t — s;y){N{v) + av){y, s) dy ds + a{t)dxu{x). 

Jo J-oo 

Defining 



(4.15) 



/oo 
-oo 



-oo 
pt r+oo 

+ ey{y,t - s){N{v) + av){y,s)dyds, 

Jo J-oo 
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following [ZH, Z4, MaZ2, MaZ3], where e is defined as in (4.8), and recalling 
the decomposition G — E+G^ + G of (4.11), we obtain the reduced equations 

/oc 
{Gu + G){x,t;y)vo{y) dy 

{^.LV) i-t Too 

/ {Gu + G)y{x,t- s;y){N{v) + av){y, s)dyds, 



-'-oo 



and, differentiating (4.15) with respect to t, and observing that ey{y, s) ^ 
as s — > 0, as the difference of approaching heat kernels. 



/oo 
et{y,t)vo{y) dy 



oo 

t P+OO 



eyt{y,t - s){N{v) + av){y,s) dyds. 

-'-oo 



We emphasize that this (nonlocal in time) choice of a and the resulting re- 
duced equations are completely different from those of Section 3.1, according 
to their respective purposes. A third possible choice has been introduced in 
[GMWZl, GMWZ2] for the study of the inviscid limit problem. As discussed 
further in [Go, Z4, MaZ2, MaZ3, GMWZl, BeSZ], a may be considered in 
the present context as defining a notion of approximate shock location. 



4.4 Nonlinear damping estimate 

Proposition 4.5 ([MaZ3]). Assuming (H0)-(H3), let vq e H^, and suppose 
that for < t < T , the H"^ norm of v remains bounded by a sufficiently 
small constant, forv as in (4.1) and u a solution of (1.1). Then, for some 
constants 61^2 > 0, for all < t < T , 

(4.18) ||t^(t)||^4 < Ce-^^*i|t;(0)||^4 + C f e-'^^'-'\\v\l, + \a\''){s)ds. 

Jo 

Proof. Subtracting the equations for u{x + a{t),t) and u{x), we may write 
the perturbation equation for v alternatively as 

(4.19) Vt + i^j df(u{x) +Tv{x,t)) drv^ — v^x — di{t)dxu{x). 
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Observing that di{dxu){x) — 0{e~^^^^) is bounded in norm for j < 4, we 
take the inner product in x of X]j=o ^x"*^ against (4.19), integrate by parts 
and rearrange the resulting terms to arrive at the inequality 

dtMUt) < -d\\dlv\\l, + C {Ml, + \a{t)\') , 

^ > 0, for C > sufficiently large, so long as ||v||i/4 remains bounded. Using 
the Sobolev interpolation 

\\vrH.<C-'\\dlv\\l+C\\v\\l. 

for C > sufficiently large, we obtain 

dt\\v\\l4t) < -Hv\\h^ + C {\\v\\l2 + \ait)\'') , 
from which (4.18) follows by Gronwall's inequahty. □ 

4.5 Proof of nonlinear stability 

Decompose now the nonlinear perturbation v as 

(4.20) v{x,t) = w{x,t) + z{x,t), 
where 

(4.21) w := UcsV, z := U^v. 

Applying U^s to (4.16) and recalling commutator relation (4.5), we obtain 
an equation 



/oo 
G{x,t;y)wo{y) dy 

[4.22) 



oo 

t poo 



Gy{x,t - s;y)Ilcs{N{v) + av){y,s)dyds 

for the flow along the center stable manifold, parametrized by w e T^cs- 

Lemma 4.6. Assuming (HO)-(Hl), forv lying initially on the center stable 
manifold A4cs, 

(4.23) \z\wr,P < C\w\jj2 

for some C > 0, for all 1 < p < oo and < r < A, so long as \w\h^ remains 
sufficiently small. 
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Proof. By tangency of the center stable manifold to S^s, we have immedi- 
ately \z\fj2 < C\w\'jj2, whence (4.23) follows by equivalence of norms for 
finite-dimensional vector spaces, applied to the p-dimensional subspace E„. 
(Alternatively, we may see this by direct computation using the explicit de- 
scription of UuV afforded by Lemma 4.1.) □ 

Proof of Theorem 1.2. Recalling by Theorem 1.1 that solutions remaining for 
all time in a sufficiently small radius neighborhood M of the set of translates 
of u lie in the center stable manifold M-csi we obtain trivially that solutions 
not originating in M.cs must exit M in finite time, verifying the final assertion 
of orbital instability with respect to perturbations not in M.cs- 

Consider now a solution v G M-csi or, equivalently, a solution w G T^cs of 
(4.22) with z = ^cs{w) G S^,. Define 

(4.24) — sup {\w\H2{l + s)-i + \w\L^ + \a{s)\{l + s)^. 

0<s<t ^ ' 

We shall establish: 

Claim. For alH > for which a solution exists with C uniformly bounded 
by some fixed, sufficiently small constant, there holds 

(4.25) CW<C2(^o + CW') for := |^;o|Lin//2. 



From this result, provided < 1/4C|, we have that ({t) < 2C2E0 
implies ({t) < 2C2EQ, and so we may conclude by continuous induction that 

(4.26) C(*) < '2C2E0 

for all t > 0, whence we obtain the stated bounds by definition (4.24). (By 
Lemma 4.5 and standard short-time existence theory, v E exists and ( 
remains continuous so long as ( remains bounded by some uniform constant, 
hence (4.26) is an open condition.) Thus, it remains only to establish the 

claim above. 

Proof of Claim. We must show that u{9 + ipi + 'ip2)^^ and |q;(s)|(1 + s) 
are each bounded by C{Eo + C(^)^); ior some C > 0, all < s < t, so long 
as C remains sufficiently small. 

By Lemma 4.6, jwolLinifa < l^'olLin^a + l^oUinifz < boUinifz + C\wo\h2, 
whence 

|wo|Lini?2 < CEq. 
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Likewise, by Lemma 4.6, (4.24), (4.4), and Lemma 4.1, for < s < 
(4.27) \U,,{N{v)+av)iy,s)\L. < CC(t)'(l + 5)-t. 

Combining the latter bounds with representations (4.22)-(4.17), taking 
g = 2 in (4.27), and applying Corollary 4.3, we obtain 
(4.28) 

POO 

\w{x,t)\LP < / G{x,t;y)wo{y)dy 



LP 



+ 



t POO 



Gy{x,t- s]y)Ilcs{N{v) + av){y,s)dyds 

^-oo 



LP 



< Eo{l + t)~^^'-p 

pt POO 

+ CC{ty / {t-sy-^^^{l + s)-'Uyds 

Jo J-oo 

and, similarly, using Holder's inequality and applying Corollary 4.4, 
(4.29) 

POD 

Ht)\ < / \etiy,mvo{y)\dy 



pt p-\-oc 

+ / / \eyt{y,t - s)\\Ilc^{N{v) + dv){y,s)\dyds 

Jo J-oo 

< |et|Loo|tio|Li + CC(i)^ / [eytlL^it - s)\flcs{N{v) + av)\L2{s)ds 

Jo 

<Eo{l + t)-'^+CC{tf [ {t-s)-^l + s)-Us 

Jo 

<c{Eo+cm{i+t)-'^. 

Applying Lemma 4.5 and using (4.28) and (4.29), we obtain, finally, 

(4.30) Iwln^it) < CiEo + C(t)')(l + t)-i 

Combining (4.28), (4.29), and (4.30), wc obtain (4.25) as claimed, completing 
the proof of the Claim and the theorem. □ 

Remark 4.7. Wc point out that the finite-dimensional part z of v is in fact 
controlled pointwise by its norm and thus by \w\h2, satisfying |^(x, t)| < 
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Ce~^l^l|2;(-, t)\L2(^x) < Ce~^l^l|u'(-, ^)|^2(2.), making possible a pointwise version 
of the argument above. This is a key point in treating the nonclassical over- or 
undercompressive cases, which appear to require pointwise bounds [HZ, RZ]. 

Acknowledgement. Thanks to Milena Stanislavova and Charles Li for 
two interesting discussions that inspired this work, and to Milena Stanislavova 
for pointing out the reference in [GJLS]. 



A Proofs of miscellaneous lemmas 

We include for completeness the proofs of earher cited lemmas that were not 
proved in the main body of the text. 

Proof of Lemma 2.2. By Taylor's Theorem, 

9{f2{t)) - gUiit)) = dg{f,mf2{t) - + omt) - A(t)i), 

as \f2{t) — fi{t)\ — >• 0. By the assumed uniform boundedness of \dg\, we 
readily obtain that / —>■ dgf is a bounded linear operator from Br, ^ Bjj. On 
the other hand \f2-f1\it) < e-''(*-*°)||/2-/i||^ for t > to implies I/2-/1I ^ 
uniformly on t > to as II/2 — /iHr? 0, and so we have 

\\9{f2{t)) - g{fi{t))-dg{f,mf,{t) - f.mu = 

^^■'^^ supe^(*-*°)o(|/2(t) - f\m = 0(11/2 - /ill,, 

t>to 

yielding the result for k = 1. The general result then follows by induction on 
k, applying the result for A; = 1 to successively higher derivatives of g. □ 

Proof of Lemma 2.4- (i) Triangulating, we have 
(A.2) 

\y{x2) - y{xi)\ = \T{x2,y{x2)) - T{xi,y{xi)\ 

< \T{x2,y{x2)) - T{x2,y{xi)\ + \T{x2,y{xi)) - T{xi,y{xi)\ 

< 9\y{x2) - y{xi)\ + L\x2 - xi\, 

where < ^ < 1 and < L are contraction and Lipschitz coefficients, yielding 
after rearrangement \y{x2) — y{xi)\ < jz:g\x2 — a;i|. 

(11) Applying Taylor's Theorem, and using the result of (i), we have 
(A.3) 

y{x2) - y{xi) = T^{x2 - xi) + Ty{y2 - yi) + o(|,T2 - xi \ + \y{x2) - y{xi)\) 
= Tx{x2 - xi) + Ty{y2 - yi) + o(|a;2 - xi\), 
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where all derivatives are evaluated at {xi,y{xi)). Noting that the operator 
norm \Ty\ is bounded by contraction coefficient < ^ < 1, we have by 
Neumann series expansion that (Id— T^) is invertible with uniformly bounded 
inverse |(Id — Ty)^^\ < (1 — 9)~^. Thus, rearranging, we have 

y{x2) - y{xi) = (Id - Ty)~^T^{x2 - xi) + o{\x2 - xi\), 

yielding the result for A; = 1 by definition of (Prechet) derivative. The results 
for A; > 1 then follow by induction upon differentiation of (2.6). □ 

Proof of Lemma 2. 7. The Lipshitz constant is bounded by max I^^A^"^ | , where 

^\{\w\le)fl{\w\le){N{t,w)l\w\)^p{\w\le)d^N{t,w)\ 
<2e\ max r p ( r ) max — h max — — 

\ hl<2£o \wY |«'|<2eo \w\ 

< 2£( max |r|p'(|r|) max +1) max |ft^,A^(i, 

the final inequality following by A^(t, 0) = 0, dy:,N{t, 0) = and the Integral 
Mean Value Theorem, or first-order Taylor remainder formula. □ 

Proof of Corollary 3.5. By Lemma 3.4, it is sufficient to prove the corre- 
sponding bounds for the purely differential operator L. By sectoriality of L, 
we have the inverse Laplace transform representations 



(A.4) 



He, := / e"*(A-L)-iciA, 



where Vcs denotes a sectorial contour bounding the center and stable spec- 
trum to the right [Pa], which by (Al) may be taken so that KF^ < uj, 
and r„ denotes a closed curve enclosing the unstable spectrum of L, with 
3f?r, > /? > 0. 

Applying the resolvent formula L(A — L)~^ = A(A — L)~^ — Id, we obtain 
in the standard way 



e*^Tlj := ! Xe^\X-L)-U\, 
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from which we obtain immediately the second stated bound, and, by a scahng 
argument [Pa], the bound 

(A.5) ||e*^n,J^,i^;,3 < ||Le*^n,||^i^^i < C(l + ri)e'^*. 

Rccalhng the standard bound ||e*^ncs||//i^//i < Ce^^, and interpolating be- 
tween I • liji and I • |ii-3, we obtain the first stated bound. □ 

Proof of Lemma 3.6. The norm in is a quadratic form, hence the map 

is smooth, and is as regular as Qq. Now 

< 3\vi-V2\h2[ sup h sup \dgo{v)\Hi], 

and sup|^|^2<5 \QQ{v)\m = 0{6^), sup^^^^^^s \dgo{v)\m = 0{6). □ 
Proof of Corollary 4-4- For definiteness, take y < 0. Then, (4.8) gives 

ey{y, t) := J2 K-^^k' {K{y + a^t, t) - K{y - a^t, t)) 

et{y, t) := ^ {{K + Ky){y + a'^t, t) - {K + Ky){y - a^k\ t)) , 

etyiy,t) := K-K' {{Ky + Kyy){y + a-^t,t) - {Ky + Kyy){y - a^t,t)) , 
where 
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denotes the standard heat kernel. The pointwise bounds (4.4)-(4.4) follow 
immediately for t > 1 by properties of the heat kernel, in turn yielding (4.4)- 
(4.4) in this case. The bounds for small time t < 1 follow from estimates 



\Ky{y + at,t,P) - Ky{y - at,t,P)\ 



y-at 



Kyy{Z, t, P) dZ 



y+at 



py-at 2 2 
Jy-\ 



I y+at 



and, similarly, 

\Kyy{y + at,t,(3) - Kyy{y - a,t,(3)\ 



at 



Kyyy{Z,t,P) dZ 



at 



< cr 



y-at 



eMt dz < Ct e Mt_ 



y+at 



The bounds for \ey\ are again immediate. 



□ 
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